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Symmetry energy for asymmetric nuclear matter at sub-saturation densities was investigated in 
the framework of an isospin-dependent quantum molecular dynamics model. Single ratio of neutrons 
and protons is compared, for the first time, with the experimental data of Famiano et al. We have 
also performed the comparison for double ratio with experimental as well as different theoretical 
resuhs of BUU97, IBUU04, BNV and ImQMD models. It is found that the double ratio predicts 
the softness of symmetry energy, which is little underestimated in single ratio. Further, the study of 
single ratio is extended for different kind of fragments, while, double ratio is for different neutron-rich 
isotopes of Sn. 
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I. INTRODUCTION 

One of the most important challenges in heavy-ion 
physics is the determination of the isospin dependence of 
nuclear equation of state (NEOS), which plays very im- 
portant role at low energy phenomena like nuclear struc- 
ture, nuclear astrophysics [l[ , fusion, cluster radioactivity 
etc; intermediate energy phenomena like multifrag- 
mentation, stopping, flow j3l426l| etc.; and at last high 
energy phenomena like pion and kaon production [27l . [2a | 
etc. The symmetry energy is found to be the prominent 
candidate to study the isospin dependence of NEOS. In 
the past years, many studies are performed on the den- 
sity dependence of symmetry energy at sub-saturation 
densities by using isotopic scaling THlll. isobaric ratio 
, single and double ratios [Ml Il9ll20| . isospin diffu- 
sion d, [2l| , isospin distillation/fractionation 0, [22| , and 
isospin migration/drift [1, 0, [1^ etc. Apart from these, 
transverse and elliptic flow of neutrons and protons are 
also considered as a good candidate to emphasize on the 
importance of density dependence of symmetry energy 
[2J-l2g. Even with the help of these studies, the exact 
determination of symmetry energy is still under the way. 

In present work, we only want to address the effect of 
symmetry energy on kinetic energy spectra of nucleons 
as well as the neutrons to protons ratio parameters. The 
later one was considered to be the ever first prominent 
candidate to extract the density dependence of symmetry 
energy. 

Before we discuss the contents of the present work, let 
us have some highlights of the single and double ratio 
in heavy ion collisions. The single ratio study in heavy- 
ion collisions have already been done by different exper- 
imental and theoretical groups [HHH^. In the ex- 
periments, near Fermi energy, Hilscher et al. [T^ found 
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that single ratio of pre-equilibrium nucleons is consis- 
tently higher than that of projectile-target system and 
it can not be explained by the Coulomb effects alone. 
Another experimental observation is the ratio of free 
neutrons and protons from two isotopic systems at 26 
MeV/nucleon. A lot of interestin g ob servations are made 
from the data. Schroder et al. [23| also systematically 
studied the spectra of pre-equilibrium neutrons and pro- 
tons in both isospin symmetric and asymmetric systems. 
Recently, at NSCL/MSU Famiano et al. [l^ measured 
the single and double ratios of free neutrons to protons 
for "2Sn 4- "^g^ and ^^^Sn + i24Sn at 50 MeV/nucleon. 
The results of double ratio of the above data have also 
been reproduced by different theoretical models, such as 
BUU97 [H, IBUU04 0, BNV [13 and ImQMD 
Even, there are a lot of uncertainties in the determina- 
tion of symmetry energy in term of different parameters 
like cross-section, symmetry energy coefficient, impact 
parameter and method of clusterization etc. 

However, no study exists in the literature where the 
comparison of single ratio of neutrons to protons is per- 
formed with the experimental data. One step ahead, the 
single ratio for the fragments is still poorly known in 
the literature. A few studies existed from the BNV and 
IBUU04 calculations for the single ratio using the inter- 
mediate mass fragments (IMF's), which is only limited 
for the small range of the kinetic energy however, 
isospin distillation/fractionation is studied up to higher 
kinetic energy by Li et al. [30j . In extension of single 
ratio to double ratio, no one has tried to compare the 
double ratio findings for experiments and theories at one 
place to see which one is the most appropriate model 
and symmetry energy form. It is also absent from the 
literature that what is the affect on the double ratio if 
we consider the series of isotopes with different isospin 
contents? From all these gaps, it seems interesting to 
perform study on the single and double ratios simultane- 
ously. 

In this paper, we focus on the comparative study of 
single and double ratios of neutrons to protons with the 
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experimental data of the MSU/NSCL collaborations [131 . 
Moreover, in addition to the comparison with the exper- 
imental data, our results of IQMD model (initially devel- 
oped by Hartnack et al. [3l|) are also compared with 
other studies of BUU97, IBUU04, BNV and ImQMD 
models. The study of single ratio is extended for differ- 
ent kinds of fragments up to higher kinetic energy, while, 
the double ratio is investigated for different neutron-rich 
systems having different isospin content. 

The article is organized as follow: we discuss the model 
briefly in Sec. II. Our results and discussions are given 
in Sec. Ill and we summarize the results in Sec. IV. 



II. FORMALISM: ISOSPIN DEPENDENT 
QUANTUM MOLECULAR DYNAMICS MODEL 
(IQMD) 

In IQMD model 0, 

US im , nucleons are represented by 
the wave packets, just like the QMD model [s^. These 
wave packets of the target and projectile interact by the 
full Skyrme potential energy, which is represented by U 
and is given as: 



U = Up + Ucoui- 



(1) 



Here Ucoui is the Coulomb energy, and Up is originated 
from the density dependence of the nucleon optical po- 
tential and is given as: 
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> = Pn + Pp, Pn and Pp are the 

neutron and proton densities, respectively. The densities 
p, p„ and Pp has the dimensions of fm~^. 

First two of the three parameters of the Eq. [5] (a and 
/?) are determined by demanding that at normal nuclear 
matter densities, the binding energy should be equal to 
16 MeV and the total energy should have minimum at 
Pq. The third parameter 7 is usually treated as a free pa- 
rameter. Its value is given in term of the compressibility: 



(3) 



The different values of compressibility give rise to a soft 
and a hard equation of state. The soft equation of state 
is employed in the present study with the parameters 
a = -356 MeV, /3 = 303 McV and 7 = 7/6 corresponding 
to isoscalar compressibility oi k = 200 MeV. In the cal- 
culations, we use the isospin dependent in-medium cross 
section in the collision term and the Pauli blocking effects 
just like QMD model [HI. The third term in the Eq.[2]is 
the symmetry potential energy for a finite nuclear mat- 
ter. The symmetry energy per nucleon employed in the 
simulation is the sum of the kinetic and potential term. 



So, the total symmetry energy is given as: 
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where, Cg^k = 25 McV from the Fermi Dirac distribution, 
which is well explained in Ref. (ssj . is known as symme- 



try kinetic energy coefficient, while, Cg 



35.19 MeV 



is parametrized on the basis of the experimental value 
of the symmetry energy, is known as symmetry potential 
energy coefficient. On the basis of 7^ value, symmetry 
energy is divided into two types with 7^ = 0.5 and 
7i = 1.5, corresponds to the soft and stiff symmetry 
energies, respectively. Note that the 7 used in the deter- 
mination of equation of state and 7^ used in the determi- 
nation of symmetry energy are the different parameters. 
The interesting feature of symmetry energy is that, its 
value increases with decreasing 7j at sub-saturation den- 
sities, while, opposite is true at supra-saturation densi- 
ties. In other words, the soft symmetry energy is more 
pronounced at sub-saturation densities, while the stiff 
symmetry energy at supra-saturation densities. 

The cluster yields are calculated by means of the coa- 
lescence model, in which particles with relative momen- 
tum smaller than Ppermi and relative distance smaller 
than i?o are coalesced into a cluster. The value of Rq and 
Ppermi for the present work are 3.5 fm and 268 MeV/c, 
respectively. 



III. RESULTS AND DISCUSSION 

In the present study, we simulate thousands of events 
for the isotopes of Sn, namely -^^^Sn -I- ^-^^Sn, -^^^Sn 
-I- i24sn and ^^^Sn + ^^^Sn at incident energy of 50 
MeV/nucleon by using the soft and stiff symmetry en- 
ergy having 7^ = 0.5 and 1.5, respectively. The colli- 
sion geometry for the study is from scmi-ccntral to semi- 
peripheral one by keeping in mind the importance of im- 
pact parameter of NSCL/ MSU collaboration's experi- 
mental results. As discussed earlier, the soft equation 
of state with an isospin dependent NN cross sections of 

CTmed = (^^ — ^''^'po) f^'^<^ employed. The single and 
double ratio is considered as a point of importance in the 
present study. The neutrons to protons ratio is among 
the first observables that was proposed a possible sensi- 
tive probe for symmetry energy prediction [Tsl . [T^ . [20| . 
This ratio is studied for the free nucleons, light charged 
particles (LCP's) (having charge number of 1 and 2) and 
intermediate mass fragments (IMF's) (having charge be- 
tween 3 and Ztot/6), where Ztot is the total charge of the 
projectile and target under study. The single ratio is just 
the ratio of neutrons to protons and is represented in the 
study by the Rn/z, while double ratio is the ratio of the 
single ratios of any two isotopes of the Sn. In order to 
study the systematics of the isospin effect, the single ra- 
tio of the isotope with more number of neutrons is always 
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mentioned in the numerator, when double ratio is calcu- 
lated. Mathematically, the double ratio is represented by 
DRn/z and is given as : 

j^neutron rich 

DRn/Z = r. (5) 

T>rieutron weak ^ ' 

A. Kinetic energy spectra 

In order to go in detail in the results from all above ra- 
tios, let us understand the kinetic energy spectra of pro- 
tons and neutrons for all type of fragments in the center 
of mass frame. The spectra of free protons and neu- 
trons are very important experimental observables that 
can provide useful information about the particle pro- 
duction mechanism and reaction dynamics. 

Fig. [T] displays the kinetic energy spectra for protons 
and neutrons at incident energy E = 50 McV/nucleon for 
semi-central geometry, while. Fig. [2]is at semi-peripheral 
geometry. The results are for neutron-rich ^■^^Sn + ^■^^Sn 
and neutron-weak system ^^^Sn -I- ^^^Sn by using the 
soft and stiff symmetry energy, respectively. The left and 
right panels in both figures are with the soft and stiff sym- 
metry energy, while top, middle and bottom panels are 
for the free nucleons, and bound nucleons inside LCP's 
and IMF's, respectively. 

It is clear from the figure that the production of neu- 
trons is more favorable for neutron-rich system [s^ ]. It 
is also true for all types of fragments as well as for the 
soft and stiff symmetry energy. This is due to the reason 
that in more neutron-rich system, the symmetry energy 
is more repulsive (attractive) for neutrons (protons) and 
hence more neutrons can be produced. On the second 
hand, the difference between yield or content of neutrons 
and protons decreases with the increasing of the kinetic 
energy. This is due to the Coulomb repulsion, which 
shifts the protons from low to high kinetic energy. The 
behavior is the same for all types of fragments. Interest- 
ingly, more neutrons can be produced with the soft sym- 
metry energy for free particles as compared to the stiff 
one. The opposite is true for the LCP's and IMF's up to 
a certain kinetic energy and after that the same trend is 
observed just like for the free particles. It is an interesting 
phenomenon and unfortunately no one has noticed this 
one. This is due to the reason that the Coulomb effects 
are stronger inside LCP's and IMF's as compared to free 
particles. However, at a sufficient high kinetic energy, 
symmetry energy dominates over the Coulomb interac- 
tions and the behavior becomes just similar to that of free 
nucleons. Although, this intersection between soft and 
stiff symmetry energy for the fragments is not so clearly 
observed from here, so, we have extended the study with 
single ratio Rn/z in the Fig. [31 

Let us move to the Fig. [2l which is displayed at semi- 
peripheral geometry. Almost, the same spectra are ob- 
served at semi-peripheral geometry, except for some ex- 
ceptions. Once again, interesting point is that yield of 



free neutrons at high kinetic energy for semi-central ge- 
ometry (Fig. [1} is higher in comparison with the semi- 
peripheral one (Fig. [2|). As we already knew that the 
symmetry energy/potential has two important functions: 
firstly, it tends to unbound more neutrons and secondly, 
it makes the neutrons more energetic than protons. Due 
to this, most of the finally observed neutrons are un- 
bounded in the very early stage of the reaction as a result 
of the nucleon-nucleon (NN) collisions at semi-central col- 
lisions. Now, symmetry energy at pre-equilibrium time is 
just shifting the more neutrons towards the high kinetic 
energy. On the other hand, at semi-peripheral geometry, 
the emission of neutrons also depend on the symmetry 
potential/energy due to the relative lack of the NN col- 
lisions. The symmetry energy makes the neutrons un- 
bound, but at relatively low kinetic energy. That is why, 
isospin effects are more pronounced at low kinetic energy 
for peripheral collisions and at high kinetic energy for 
central collisions. These results are also consistent with 
those in Ref. 0. When one moves from LCP's to IMF's, 
the content of neutrons for neutron-rich system is lower 
than the neutron-weak system at high kinetic energies 
for semi-central as well as semi-peripheral geometries. 
This is due to the reason that at high kinetic energy the 
yield of free nucleons is higher as compared to fragments 
and it will result in more production of free neutrons for 
neutron-rich system as compared to the fragments. 

B. Single ratio 

In order to make sure about the above discussion from 
Fig. [T] and Fig. [5J it is interesting to investigate the 
single ratio (Rn/z) of neutrons to protons for free nu- 
cleons, LCP's and IMF's, which is shown in Fig. [3] for 
neutron-rich system ^^^Sn + ^"^^Sn and neutron- weak sys- 
tem ^^^Sn -I- -'^^^Sn by using the soft and stiff symmetry 
energy. The left and right panels are at semi-central and 
semi-peripheral geometries, respectively. As is expected 
from Fig. [1] and Fig. [21 Fig. [31 depicts the results as 
follow: 

• The isospin effects for more neutron-rich system are 
stronger and it is consistent with Ref. Q and with 

Figs, m n 

• Rn/z decreases with the kinetic energy for all 
types of fragments at semi-central as well as semi- 
peripheral geometries. 

• For free nucleons. the isospin effects are stronger 
at high kinetic energy for semi-central geometry, 
while, the same is true at low kinetic energy for 
semi-peripheral geometries. It is also explained ear- 
lier in Ref. @ . 

• The increase in the neutrons to protons ratio for 
neutron-rich system at sufficient high kinetic en- 
ergy is due to the repulsive nature of the symmetry 
energy for neutrons. 
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FIG. 1: (Color online) The kinetic energy spectra in the center 
of mass system for neutrons (solid symbols) and protons (open 
symbols) from the free nucleons (a,b), LCP's (c,d), and IMF's 
(e,f) at semi-central geometry ( b = 2 fm ) of ^'^^Sn + ^^^Sn 
(blue circles) and "^Sn + "^Sn (black squares) collisions at 
E = 50 MeV/nucleon by using the soft (left) and stiff (right) 
symmetry energy, respectively. 



• Let us discuss the single ratio for fragments. R^/z 
of IMF's is earlier studied by the Catania group 
using the BNV fisf and Texas group using the 
IBUU04 model . Both models have different ap- 
proaches for the symmetry energy and hence the 
results are little different from each other. In the 
BNV results, ratio decreases at low fragment ki- 
netic energy and then increases at high kinetic en- 
ergy for neutron-rich system with the stiff symme- 
try energy. On the other hand, in the IBUU04 cal- 
culations, the ratio is found to decrease with frag- 
ment kinetic energy. However, both the calcula- 
tions have the same behavior with the soft and stiff 
symmetry energy. But, both groups have limited 
their study only to the relative low kinetic energy 
and were not able to investigate the cross-over phe- 
nomenon of symmetry energy, which takes place at 
higher kinetic energies and discussed in detail in 
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FIG. 2: (Color online) Same as in Fig. [U but for the semi- 
peripheral geometry (b= 6 fm). 



this study. 

• The large isospin effects are observed with the soft 
symmetry energy for free nucleons along whole 
range of the kinetic energy [23|, while cross-over 
happens for the LCP's and IMF's at certain ki- 
netic energy. Below the cross-over kinetic energy, 
the stiff symmetry energy produces larger neutrons 
to protons ratio and after the cross-over, it is true 
with the soft symmetry energy and behaves like 
just for free nucleons. Recently, Harmann et al. 
p35| displayed the data for single ratio from the 
IMF's below the cross-over kinetic energy. This 
data (not shown here) is favoring the soft symme- 
try energy in our studies with the IQMD, however, 
with the BNV, there data is well explained by the 
stiff symmetry energy. If one sees them carefully, 
we can find that the soft symmetry energy is more 
soft and stiff symmetry energy is less stiff in the 
BNV as compared to the IQMD and ImQMD mod- 
els. In other words, the stiff symmetry energy from 
the BNV and the soft symmetry energy from the 
IQMD /ImQMD lies between the stiff symmetry en- 
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FIG. 3: (Color online) The ratio of neutrons to protons (a,b) 
at freeze-out time from free nucleons, LCP's (c,d) and IMF's 
(e,f) as a function of kinetic energy at semi-central (left) and 
semi-peripheral (right) geometries by using the soft (solid 
symbols) and stiff symmetry energy (open symbols). The 
vertical lines in the plots of LCP's and IMF's represent the 
kinetic energy at the cross-over points of the soft and stiff 
symmetry energy. Blue circle represents for ^'^'^Sn -I- ^^^Sn 
and black squares for "^Sn + "^Sn at E = 50 MeV/nucleon. 



ergy from the IQMD/ImQMD and the soft sym- 
metry energy from the BNV calculations. It means 
that the data is favoring almost the same symmetry 
energy from both the models. 

• The cross-over kinetic energy is at higher value for 
more neutron-rich system and increases with the 
size of the fragments (i.e. from LCP's to IMF's). 

• The cross-over value of the kinetic energy also 
raises when one moves from semi-central to semi- 
peripheral geometries. This value is more affected 
for the more neutron-rich system. 

• It is also clear from here that gas phase (free nucle- 
ons) is significantly enriched in neutrons relative to 
the liquid phase or fragments that are represented 
by the bounding nuclei. The phenomenon is known 
as isospin distillation/fractionation and is discussed 
many times in the literature only in term of the free 



FIG. 4: (Color online) The comparison of neutrons to pro- 
tons ratio from free nucleons, for the systems ^^'^Sn -I- ^^^Sn 
(top) and ^^''Sn + ^^^Sn (bottom) at E = 50 MeV/nucleon 
and impact parameter b < 5, with the experimental data of 
MSU/NSCL collaborations The filled and open circles 

represent the soft and stiff symmetry energies, respectively. 



and bound nucleons [6[. More interesting results 
are expected for isospin distillation if one tries to 
study in term of different kind of fragments. 

The theoretical results become more interesting and 
useful if one compares the results with experimental data. 
In the Fig. |4l we have, for the first time, compared the 
results of single ratio of neutrons to protons from free 
nucleons for neutron-weak system ^^^Sn -|- -'^^^Sn (in top 
panel) and for neutron-rich system ^^^Sn -f- ^^^Sn (in bot- 
tom panel) at 50 MeV/nuclcon with the experimental 
data [1^. The behavior of single ratio results for free 
nucleons is explained in the Fig. [31 The conclusion from 
the figure is that (1) more Rn/z is observed for more 
neutron-rich system, which is also predicted by theoreti- 
cal predictions, and (2) R^/z shows increment at higher 
kinetic energy, especially for ^^"^Sn -t- ^^'^Sn. It indicates 
that the theoretical results are consistent with the exper- 
imental one. 

The results are in good agreement with the soft sym- 
metry energy except at very low and very high kinetic 
energy. The difference between soft and stiff symme- 
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try energy results for the neutron- weak system is almost 
comparable to the error bar, while, for neutron-rich sys- 
tem, the difference has a great importance over the error 
bar. In other words, the error bar of the theoretical re- 
sults with the soft symmetry energy covers the error bar 
of the experimental data for both the systems under con- 
sideration. The difference at high kinetic energy between 
theoretical and experimental results for neutron-rich sys- 
tem is due to the large uncertainty in the measurement of 
the Rn/z- By using the single ratio observable, one can 
reach at a partial conclusion that the asymmetric nuclear 
matter favors the soft symmetry energy at sub-saturation 
densities, which is also consistent with the other findings 
in the literature 

As we have observed, the single ratio mixes the symme- 
try energy with Coulomb effects throughout the kinetic 
energy range. In order to minimize the Coulomb effects 
and systematical error, it is reasonable to study the dou- 
ble neutrons to protons ratio for the isotopes of the same 
element. This is also studied in the literatures with only 
two isotopes [l^ [H, [l^ [131 ■ No one has tried to inves- 
tigate the effect of double ratio on a series of isotopes in 
the asymmetric nuclear matter so far. 
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C. Double ratio 

In the present study, we consider reactions between 
three isotopes of Sn and observe the relative effect of 
these isotopes on the double ratio and symmetry energy. 
The pairs arc as follow: ^^^Sn + ^^^Sn and ^^"^Sn + ^^^Sn, 
i24Sn + i24sn and "^Sn + "^gn^ 132^^ + 1323^^ 

^^^Sn -f ^^^Sn. The three pairs having the difference of 
8,12 and 20 neutrons, respectively. The universal behav- 
ior for double ratio is observed with the kinetic energy, 
i.e. with the increasing of kinetic energy, the double ra- 
tio is found to increase for all the three sets of isotopes 
which we have plotted in the Fig. [5] The increase in 
the double ratio is due to the effect that now energetic 
nucleons are more affected by the symmetry potential, 
which are already suppressed by the Coulomb repulsion 
in the single ratio results. The effect of symmetry en- 
ergy on the double ratio is just like the single ratio, i.e. 
larger value with the soft symmetry energy as compared 
to the stiff one. Moreover, double ratio goes on increas- 
ing with the increase of the neutron difference between 
the pairs discussed above, or in other words, the dou- 
ble ratio from free nucleons goes on increasing with the 
initial-state double ratio of the systems from three differ- 
ent pairs of isotopes of Sn. 

This increase is due to the effect that the more the 
number of neutrons, the more repulsive the symmetry 
energy for them. The Coulomb effects are already can- 
celled by taking the double ratio. Hence, the results are 
just like that as expected. The double ratio is found to 
be weakly sensitive towards the collision geometry. How- 
ever, a little increase is observed at semi-peripheral geom- 
etry compared to semi-central one at high kinetic energy. 



FIG. 5: (Color online) Free neutrons to protons double ra- 
tio, at semi-central (left) and semi-peripheral (right) geome- 
try with the soft (top) and stiff (bottom) symmetry energy, 
as a function of kinetic energy at the incident energy E = 50 
MeV/nucleon. The different lines in the figure are the double 
ratio from different pairs: Solid line for ^^^^Sn + ^^^Sn and 
^^■'Sn + i^^'Sn; Dashed line for ^^^Sn + ^'^'^^n and ^^^Sn + 
^i^Sn; Dash-dot line for ^^^Sn + ^^^Sn and "^Sn + "^Sn. 



This is true with the stiff as well as the soft symmetry 
energy. 

The double ratio has been studied many times in last 
couple of years by different groups with the help of the 
BUU97, IBUU04, BNV and ImQMD models and com- 
pared with the experimental results. Even so, we are still 
far away from the exact conclusion about the symmetry 
energy form. We have, along with all the possible re- 
sults in the literature, compared the double ratio with the 
IQMD model in Fig. El L et us start with very first com- 
parison of the BUU97 [T^l- The results were very close to 
the experimental one, but, the reaction conditions were 
different. Firstly, in the BUU97 calculations, the inci- 
dent energy was 40 MeV/nucleon, not 50 MeV/nucleon, 
just like the experimental one. Secondly, data set is only 
for the transverse emission, while in the BUU97 calcula- 
tions, the nucleons used are emitted in all the directions. 
Move one step ahead to the IBUU04 results [l^, where 
the symmetry energy is introduced with the help of mo- 
mentum dependent interactions, the results are very far 
from from the experimental data. The same is true for 
the BNV calculations performed by the Catania group in 
2007 171. The most closeness between the data and the 
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IV. SUMMARY 
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FIG. 6: (Color online) The comparison of free neutrons to 
protons double ratio at E = 50 MeV/nucleon and 6 < 5 with 
the MSU/NSCL data, BUU97, IBUU04, BNV, and ImQMD 
simulations. 



calculation is observed by the ImQMD model in 2009 |19| . 
They found that the results with 7^ ~ 0.75 are best fit 
with the experimental data for impact parameter b < 2 
fm. In the present study, we have performed simulations 
for 6 < 5 fm and for the angular cuts, as mentioned in 
the experiments, with the soft and stiff symmetry energy 
and displayed the theoretical results over the whole range 
of the kinetic energy. The Fig. [6] clearly indicates that 
our results are very close to experimental data. 

If we see the comparison of theoretical and experimen- 
tal results from single (Fig. |4]) and double ratio (Fig. |6]) 
results, it seems that single ratio results require 7^ < 0.5 
to explain the data, and while the data is well explained 
by the 7^ = 0.5 for double ratio. This is due to the reason 
that single ratio suffers the effect from the Coulomb in- 
teractions in addition to symmetry energy. As our main 
purpose is to extract the symmetry energy, where double 
ratio can act as a better candidate rather than the sin- 
gle ratio. In conclusion, the results of double ratio can 
be very well explained by the soft symmetry energy with 
7; = 0.5 in comparison with single ratio, where the data 
is little underestimated by the theoretical predictions. 



In summary, we have performed a detailed analysis for 
kinetic energy spectra of free nucleons and bound nu- 
cleons inside fragments as well as the ratio parameters 
for the three reaction channels of Sn-isotopes at E = 50 
MeV/nucleon via multifragmentation. The kinetic en- 
ergy spectra of protons and neutrons from free nucleons 
and all type of fragments in the center of mass frame 
shows that the content of neutrons is more favorable for 
neutron-rich system since the symmetry energy becomes 
more repulsive (attractive) for neutrons (protons) and 
hence more neutrons can be produced. In addition, the 
difference between yields or contents of neutrons and pro- 
tons decreases with the increasing of the kinetic energy, 
which can be explained by the Coulomb repulsion shifts, 
making the protons from low to high kinetic energy. In- 
terestingly, more free neutrons can be produced with the 
soft symmetry energy as compared to the stiff one. 

From the single ratio of free neutrons to protons, it 
decreases with the kinetic energy for all types of frag- 
ments at semi-central as well as semi-peripheral geome- 
tries. However, the increase of the ratio from free nu- 
cleons for neutron-rich system is observed at sufficient 
high kinetic energy, this can be explained by the repul- 
sive nature of the symmetry energy for neutrons. For 
single ratios of LCP's and IMF's, we noticed a transi- 
tion at certain kinetic energy between the soft and stiff 
symmetry energy, while no transition for the free nucle- 
ons or gas phase. Below the cross-over kinetic energy, 
the stiff symmetry energy produces larger ratio of neu- 
trons to protons and after the cross-over, it is true with 
the soft symmetry energy and behaves like just for free 
nucleons. This transition is also found to be strongly 
dependent on the isospin of the colliding partners, size 
of the fragment and weakly dependent on the collision 
geometry. Moreover, isospin distillation is also observed 
when one moves from the gas phase to liquid phase. It is 
further more interesting to study the isospin distillation 
in term of different kind of the fragments as compared to 
consider bound fragments as a single liquid phase. 

The comparison of the theoretical results of single and 
double ratios with the experimental data emphasizes on 
the softness of the symmetry energy at sub-saturation 
densities, which is yet uncertain at the supra-saturation 
densities. Although, the single ratio study underesti- 
mates the data a little as compared to double ratio 
for the same stiffness of symmetry energy (7^ = 0.5), 
which reflects the double ratio is a relative good candi- 
date for density dependence of symmetry energy at sub- 
saturation densities because of the cancel of Coulomb 
effect between two systems. Of course, the magnitude of 
double ratio of neutrons and protons from free nucleons 
strongly depends on the initial double ratio of the sys- 
tems. It gives us an indication that it is better to study 
the isospin physics with a pair of ^■^^Sn -f ^'^^Sn and ^^^Sn 
+ "2Sn. 
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